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A HIGH PRESSURE HWULHTBD M01£CULAR BEAM HASS 


SPPCTROHETRIC SAMPLING SYSTEM 

by Carl k. Steams^ Fred J* Roh?.» 
George C. Fryburg^ and Robert A. Miller 

Levis Research Center 


SUMMARY 

The current state of understanding of free**jet high pressure sampling 
is critically reviewed and oodifications of certain theoretical and empir- 
ical considerations are presented. A high pressure free-Jet expansion, 
modulated csolecular beam, mass spectrometrlc sampling apparatus was con- 
structed and this apparatus is described in detail. Experimental studies 
have demonstrated that the apparatus can be used to sample high tempera- 
ture systems at pressures up to one«atmosphere. Condensible hi^ temper- 
ature gaseous species have been rou&nely sampled and the mass spectro- 
metric detector has provided direct identification of sampled species* 
System sensitivity is better than one tenth of a part per million* Ex- 
perimental results obtaintd with argon and nitrogen beams are presented 
and compared to theoretical predictions* These results and the respec- 
tive comparisons are taken to indicate acceptable performance of tne sam- 
pling apparatus. Results ate also given for two groups of experiments re- 
lated to hot corrosion studies* The formation of gaseous sodium sulfate 
in doped methane-oxygen flames was characterized and the oxidative vapor- 
ization of metals was studied in an atmospheric pressure flowing gas sys- 
tem to which gaseous salt partial pressures were added* 


INTRODUCTION 

The apparatus described herein was built to facilitate sampling and 
characterization of gaseous molecular species from various sources at 
pressures up to one atmosphere. The need for the apparatus derived from 
our desire to study certain chemicals reactions considered fimdamental in 
the hot corrosion process which can occur in gas turbine engines (refs. 1 
to 3). Envisioned experiments would involve the sampling of doped flames 
and flowing gas reaction tube systems* Species to be directly identified 
and measured would include stable gases, reactive intermediates associ- 
ated with flames, and low vapor pressure metal-containing high tempera- 
ture molecules* Therefore, prerequisites for the apparatus were: (1) cap- 

ability of sampling systems at atmospheric pressure, (2) capacity to sam- 
ple condensible high tenqperature gaseous species, (3) mass spectrometrlc 
detection to provide direct identification of sampled species and to meas- 
ure tneir concentration, (4) sensitivity of parts per million or better, 
and (5) ability to maintain the chemical and dynamic integrity of the 
species being sampled* 


National Research Council Resident Research Associate* 
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Because loass spectrometers are relatively low pressure devices » the 
apparatus design problem centered on extracting a sample from the high 
pressure source and transferring it to the mass spectrometer without 
altering the sample's chemical characteristics in an unpredictable way. 
Batch and microprobe techniques are precluded for samples containing 
reactive and condensible species. However, the so-called high pressure 
free-Jet expansion sampling technique appealed to satisfy our prerequi- 
sites, In this technique one samples (continuously or intermittently) 
from high pressure through an orifice into a vacuum to produce a super- 
sonic molecular beam which is compatible with and can be analyzed by a 
mass spectrometer operating in a low pressure region. This sampling 
technique has been successfully applied by a number of investigators in 
quite diverse research areas (refs, 4 to 12 and references cited therein). 

Our apparatus was designed in consultation with Drs. F, T. Greene 
and T, A. Milne of the Midwest Research Institute, The gas sample enters 
a four stag^ differentially pumped apparatus through an orifice in the 
apex of the sampling cone. The gas expands as a free-jet in the first 
vacuum stage and becomes a supersonic "molecular beam," The central core 
of the t>eam is passed by a conical skimmer into the second vacuum stage 
where the beam is modulated by a rotating segmented disk. After passing 
through another differential pumping stage the modulated beam enters the 
fourth very* low pressure vacuum stage and the ion source of a quadrupole 
mass spectrometer. 

Detailed evaluation of the apparatus we constructed established that 
it satisfied the prerequisites and design criteria. The system has been 
successfully used to investigate (1) the formation of sodium sulfate in 
doped flames and flowing gases at atmospheric pressure (ref. 13), and 
(2) volatile species emanation from metal samples exposed at elevated tem- 
peratures to one atmosphere pressure of oxygen containing partial pressures 
of NaCl(g) and H20(g) (refs. 14 and 15), The sampling apparatus has poten- 
tial to be applied to a variety of problems where it is desirable to sam- 
ple systems at pressure up to one atmosphere. 


THEORETICAL ASPECTS OF HIGH PRESSURE SAMPLING 
Beam Forma^on 

When san^ling a gas system through an orifice that is not in thermo- 
dynamic equilibrium with the source gas, the orifice size should be so 
selected to insure that the sample represents the bulk source gas and not 
a stagnant boundary layer near the orifice. For orifice Knudsen numbers, 
Kuq, less than about 10“^, boundary layer effects are insif.ru ticani 
(refs, 16 and 17) and therefore the orifice diameter, D , should be sized 
to satisfy the criterion KOq < lO*^ for the source gas and conditions 
of pressure, Pq, and temperature, Tq. The orifice Knudsen number is re- 
lated to the orifice diameter through the mean free path length, of 
the source gas molecules (ref, 18), i.e,. 
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( 1 ) 


VAien gas flows through an orifice into a lover pressure region^ the 
type of gas flow (effusive » transitional, slip, or continuum) is dictated 
by the value of Kn^. For KHq 10*^ the source gas expands into the 
low pressure region oy continuum flow (ref. 17). The flow produces a 
supersonic Jet if the pressure ratio. across the orifice is greater than 
about two for most gases (ref. 19). (For most beam systems the pressure 
ratio is usually made as high as practical, >10^, to minimize scattering 
by background gas in the low pressure region. ) The expansion is termed 
**free-jet" if the orifice has no flow controlling constraints on its down- 
stream side. Because a finite pressure exists in the low pressure region, 
the supersonic flow field is bounded by a shock system consisting of a 
barrel shock and a Mach disk or normal shock front. The flow field and 
shock system are shown schematically in figure 1. The size and shape of 
the shock system has been studied by Bier and Schmidt (ref. 20 and ref- 
erences cited therein) who showed that these parameters are a function of 
the pressure ratio across the orifice. Ashkanas and Sherman (ref. 21) 
have established that the location of the Mach disk is independent of the 
source gas properties. The location of the Mach disk is given by the 
equation 


% 



( 2 ) 


where is the distance downstream from the orifice to the plane of 

the Mach disk and is the pressure in the region downstream of the 

orifice. Thus we see that the higher the pressure ratio, the further 
downstream the jet proceeds before being interrupted by the normal shock. 


Continuum Flow Beam Gas History 

The gas history (time or distance-temperature-pressure-density rela- 
tions) of a free-jet expansion can be approximated as an isent topic ex- 
pansion where the gas leaves the source as a classical bulk fluid charac- 
terized by thermodynamic parameters. In the continuum flow region mole- 
cules collide with each other in such a wa> that the mass flow becomes 
highly directed and forms a Jet. In the flow process the random thermal 
energy (enthalpy) of the colliding gas molecules is converted to directed 
kinetic energy and adiabatic cooling of the gas takes place. As the ex- 
pansion proceeds downstream, the jet density and temperature decrease 
until a region is reached where collisions become so infrequent that essen- 
tially collisionless flow prevails. In this nearly colllsionless flow 
regime a supersonic molecular beam is formed. 
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A knowledge of the aerodynamics of the supersonic beam is necessary 
if one expects to estimate the time involved in the beam formation process 
and the gas history. Extensive research by numerous Investigative groups 
has established an almost complete understanding of the aerodynamics. 
Excellent reviews and detailed analyses are contained in the literature 
(refs* 19, 21 to 25 and references cited therein). Our analysis of the 
beam formation process as presented here is not intended to start with 
first principles or be all inclusive* Instead we will deal only with 
those aspects needed to provide a general background for understanding 
the general principles of the sampling system* 


The one-dimensional compressible flow analysis of a perfect gas with 
complete relaxation (equilibrium), as provided by fluid mechanics (ref. 26), 
yields the fundamental isentropic expansion equations relating the expanded 
gas properties to the flow Mach number, M, and the source gas properties 


n 




2il/(l-Y) 

M 


- ", 


(T - 1) 


y/(y-1) 


(3) 


(4) 



I 




n 





^1/(Y-1) 


(5) 

( 6 ) 


The density n, pressure P, and temperature T refer to source condi- 
tions when subscripted by zero and to expanded gas conditions at M when 
not subscripted. The specific heat ratio for the gas is denoted by y 

For incomplete relaxation, an effective y may be used. 

Mach number is defined as the ratio of the local flow veJocity, V, to the 
local speed of sound in the fluid, c, i.e. , 

V 

M = - (7) 

c 


Ashkenas and Sherman (ref, 21) have derived an expression relating 
the centerline flow Mach number to the distance X downstream in the 
expansion from the orifice. Their expression is 


M - A ^ ^0 

^ ^ - d: 


V-l 


±1\ 

- 1 ' 


Y-1 


Dq/ 


( 8 ) 
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where A and (X^/D.) are constants whose value depends on y ^ shown 
In Table I. From equation (8) it is seen that it is convenient to express 
distance downstream from the orifice in terms of orifice diameters (X/l^). 
For distances (X/D.) > ^ . the centerline Mach number can be approximated 
by 0 



(9) 


for values of y ^ 5/3, 7/5, or 9/7. For Mach nuDd>ers large in comparison 
to unity, Knuth (ref. 27} suggests the following equation for all specific 
heat ratios: 



Uhile Ashkenas and Sherman's expression (eq. (8)> has been shown 
experimentally to be quite accurate at several orifice diameters down- 
stream (refs. 21 and 28) , it has been shown (ref. 28) and is now recog- 
nized that this equation does not apply in the transonic region near the 
orifice. The transonic flow problem is still unsolved but characteriza- 
tion of this flow has been approximated by various sink flow models 
(refs. 25 and 29), by extrapolation of Ashkenas and Sherman's results 
(ref. 30), and by empirical expressions derived from e^erimental meas- 
urements (refs. 19 and 31). We have chosen to characterize the flow from 
the subsonic region upstream of the orifice to the hypersonic region down- 
stream by smoothly extrapolating equation (8) for M < 4 while assuming 
sonic flow (M » 1) at the orifice plane. Although this extrapolation is 
not based on a physical model, we believe that the results are an adequate 
approximation for our purposes. 

The extrapolation is performed by fitting the Mach number of a third 
degree polynomial in (X/D^). The polynomial 

M = Sq + aj(X/DQ) + a2(X/DQ)^ + a^fX/Dp)^ (11) 

was fit in two ranges of (X/Dq) by specifying certain conditions for each 
range. The first range was 0 <. (X/D^) ± 

the distance for M a 4 as calculated by equation (8). In this range 
the specified conditions were: (1) the flow is choked by the orifice, i.e., 
M ® 1 at (X/Dq) = 0} (2) a point of Inflection exists at the orifice, 
i.e., at (X/Dq) = 0 the second derivative, d2M/d(X/DQ)2, is zero: (3) at 
M = 4 the polynomial and equation (8) collocate , ’•>) i " A the 

pol}momial and equation (8) have the same slope. Inspection of the poly- 
nomial reveals that condition (1) sets Sq « 1.000 and condition (2) 
fixes &2 ■ 0,000, 
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The second range for fitting the polynomial extended from (X/Dq) ■ 0 
to some distance upstream from the orifice designated (X/Dq)^ which Is 

the root of the polynomial for the second range. Conditions specified 
for fitting in this range were: (1> the flow is choked by the orifice, 

(2) the polynomials for the first and second range collocate at the ori- 
fice, O) the polynomials have equal slopes at the orifice, (4) the sec- 
ond polynomial has an inflection point at the orifice, and (5) at (X/Dq)^ 
the slope is zero. 

For all conditions specified, the following general polynomials 
result for the two ranges 


Range 1: M » 1 4 s^CX/Dq) 4- a 3 (X/DQ>^ (12) 

Range 2: M » 1 4 a^CX/D^) - — Sj^X/Dq)^ (13) 

3 

The root of the second polynomial is (X/Dq) = - y, and the polynomials 

r a 

obtained for y = 5/3, 7/5, and 9/7 are presented in Table 11 . The poly- 
nomials and equation (8) were used to calculate the Mach^^umber against 
distance plots shown in figure 2 < These curves show, as ^ initially 
assumed, that the gas expansion actually begins slightly upstream of the 
orifice and this is in substantial agreement with experimental measure- 
ments and results obtained by other approximations r 

The results presented in figure 2 provide a characterization of the 
aerodynamics of the free-jet expansion which can be used to estimate the 
gas history. The hydrodynamic flow velocity is given by 

d(X/D^) 

» • ' ”o — i — 

and the speed of sound is a furjctlon of ^he local gas temperature, T, and 
molecular weight, M ^ref. 18>, 


r = 



(15) 


where R is the gas constant (8.3143 J K *-mci ^ *r 8.3i43xl0^ erg K 
mol“^) . For the isentropic expansion the local gas temperature of the jet 
can be expressed in terms of the source gas temperature (by eq (5)) giving 



( 16 ) 


or 
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c * Co(l + ^ (17) 

where Cq Is the speed of sound In the source gas. From equations (14) 
and (17) it follows that 

Integration of this equation yields the time-distance profile for the 
expanding gas* The profiles obtained by numerical integration of equa- 
tion (18) for the specific heat ratios 5/3, 7/5, and 9/7 are shown in 
figure 3* In this figure time is plotted as the dimensionless flow time 
parameter CQt/D^ and distance is scaled in terms of orifice diameters* 

Combining the data of figures 2 and 3 with the appropriate isentropic 
expansion relations (eqs* (3) to (6)) allows one to determine the tempera- 
ture, density, or pressure at any point or time in the expansion. These 
various relations are shown graphically in figures 4 and 5 where the 
parameters ratioed to source conditions are plotted against time and dis- 
tance, respectively. While these figures clearly indicate the rapidity 
of the expansion, a more tangible appreciation of the time involved can 
be obtained by considering an actual gas on a real time basis. The ex- 
pansion history of argon source gas expanding through a 0 025 centimeter 
diameter orifice is shown in figure 6, Here wf, see that major reductions 
in temperature, density, and pressure take place in only a few microsec- 
onds. This feature of the free-jet expansion accounts for the applica- 
bility of this technique to sampling (refs. 29 and 32). 


Beam Gas History Beyond Continuum Flow Region 

From the preceding discussion one might intuitively presume that at 
sufficiently large distances from the jet source the density becomes so 
low that there are essentially no collisions. Theoretical and experimental 
studies have shown that there is a limiting region in the flow field beyond 
which further cooling and acceleration are precluded by the rarefied nature 
of the flow (ref. 19). Thus there is a transition from the collision dom- 
inated continuum flow to a nearly ccllislonless molecular flow regime. Ex- 
tensive effort has been directed to understanding the transition flew 
regime and as pointed out by Anderson (ref. 19) there is a hierarchy of 
theories. This hierarchy is based on the degree of accuracy and detail 
provided by the theories. The direct simulation Monte Cctrlo method as 
applied by Bird (ref. 33) has provided an almost complete understanding of 
the transition flow regime but a detailed analysis of this regime is beyond 
the scope of our presentation. For our purposes it is only necessary to 
conq)rehend certain aspects of transition flow. Tc this end we will con- 
sider a simplified treatment which has proven to be adequate tor predicting 
some jet properties 
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Anderson and Fenn (ref, 34) measured velocity distributions of molec- 
ular beams extracted from free Jets. They found that the transition tC 
free molecular flow could for convenience be compared to phefiamenon occurs* 
ring in effusive flow. For effusive flow the velocity distribution in the 
axial direction does not change with increasing distance from the molecu- 
lar beam source. The axial (parallel) molecular velocity distribution is 
thus said to be "frozen in." Conversely, the radial (perpendicular) dis- 
tribution approaches zero with increasing distanoe from the source. (This 
is a consequence of the fact that the molecules travel in straight line 
paths and there are no collisions.) For the effusive source the distribu^ 
tion of velocities is anisotropic with a wide distribution in rhe parallel 
direction and a narrow distribution in the perpendicular direction. 

For continuum flow the distribution of velocities is isotropic with 
equal parallel and perpendicular distributions. No freezing of the par- 
allel distribution takes place in continuum flow and the velocity spread 
of the distribution decreases with increasing distance from the jet source 

In considering expanding source flow it is useful to define a par- 
allel temperature, T||, and a perpendicular temperature, These are 

the temperatures of A ‘ two-dimensional Maxwellian distribution function 
which describes the local distribution of molecular speeds in the radial 
direction and in a direction perpendicular to the radial direction. For 
continuum flow where the speed distribution is isotropic, Tj^ * Tjl. The 
transition flow problem then can be considered as one of determining how 
Tj I and Tj vary in the expansion as the transition regime is passed 
through. 


Terminal Mach Number 

In the "free molecular flow" region of a free-jet expansion Anderson 
and Fenn (ref. 34) showed that a freezing of the parallel distribution 
occurred. They found nonisotropic distributions where the perpendicular 
direction had a narrower distribution width than that of the parallel 
direction. Their model of the transition regime yielded an expression 
for a terminal Mach number, Mq«, that was in substantial agreement with 
their experimental observations and those from other laboratories. The 
terminal Mach number expression was derived by considering the point in 
the expansion beyond which the temperature cannot change as rapidly as 
required for isentroplc continuum behavior. The Mach number at this point 
M2, represents the value where collision domination of the jet flow ceases 
The terminal Mach number then represents the asymptotic value which the 
Mach number approaches at infinite distances downstream from the source. 
Anderson and Fenn*s terminal Mach number can be expressed fn the general- 
ized form 




(19) 


where A Is the constant in Ashkenas and Sherman's Mach number expression 
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ax^ the distacK:# firom the source orifice to where is YeMheel c^ he 

celcolated from equation (9) • 


£eam Intensity 

Thus far in our atialysis of the beam formation process mo tensidera- 
tion has been giren to the question of beam intensity. From a practical 
point of view, intensity is a key parameter if free*jet expansion tech- 
niques are to be used in sampling. With a conventional effusive or oven 
beam source, the maximum attainable intensity in the collimated molecular 
beam is limited by two factors for any given oven temperature. These 
factors are: (1) the effusion through the orifice and (2) geometrical 

factors involved in obtaining a collimated beam from (he random Initial 
velocities of the effusing molecules. Kantrowitz and "Orey (ref. 95) sug- 
gested that these limitations coiald be oveacone by placing a collimating 
alement in the flosj^from a sonic nozzle. Their analysis indicated that 
the mass notion from a sonic orifice provides an initial collimation which 
Impieoves the flow rate and geometrical factors so that a net beam intensi^ 
fieation results. Their work provided much of the Impetus for development 
of Che free-jet expaasion Sampling technique. 


(a) Continuum Flow at the Skinmer 

The Kantrowitz and Grey analysis is based on four critical nsaump* 
tions: (1) the continuun flow field extends to the first collimating ele- 

dent (hereafter called the skiimner) ; <2) the continuun flow is isentnopic; 
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t 


tbt £l«w at tha s^ionaar is parallel aad undisturbed by the presence 
of the skinmat* (4) the flo9 dofnatream ftom the skinoMP' is collision' 

lass. On the basis of these assumptions, and by neglecting the component 
of the random molecalak thermal velocity in the direction of the mass 
flow. Kafitvowl%z a«d Grey dbrived expressions for the velocity dlstilbu- 
Cion and "precolllmatlon gain" In beam Intensity, Parker et al. (ref. 36) 
made A similar analysis without making the mathematical slmpliflration of 
Simlect^ag tAte rasidofa thermal velocity component parallel to the beam. 

For Mach aumbefs greater than three the results of both analyses are 
essentially the same. 


V 

Tk» axial differential number density can be e|ipressed as 


dn(u) 2 

a u 

du 



(24) 


vhefe u is th# molecular velocity » V is the flow velocity » and a Is 
the most probable thermal velocity at the translational temperature T 
trtilch is given by the Isent topic flow lelation (eq. (5)). The most prob- 
able velocity of the supersonic molecular beam, is obtained by deter- 

mining the maximum for the distributi^ given by equation (24), 


u 

mp 



(25) 


Thus it is apparent that for the supersonic beam the distribution is nar- 
rower than the Mcjcwellian distribution for an effusive beam where K « 
and V « 0. Also^ the supersonic distribution is seen to move to higher 
velocities and become narrower as M increases. At the limit of infi- 
nite Mach muid>er the distribution is monoentrgetlc with velocity V, 


The theoretical equation for the centerline flux at a detector, 
located a distance i downstream from the skimmer, derived by Parker et 
al. is given by 




M_(3 + vM^Hl + erf/vM^/2> ^ 
s s s 


(2 + /M )exp 
s 


( 4 ) 


(1 + (y - DM^/ 21 ^''^ 


(26) 


vhata ttg is the density at the skimmer, Ag is the area of the skimmer 
orifice, and Mg is the Mach number at the skimmer This equation is 
baaed on the four idealization assumptions of Kan^rowitz and Grey. For 
Mg > 3, the centerline flux at the detertor ^an be approximated by 



M (vM^ + 3) 
& s 


1 ^ (y - DMV? 

s 




1/2 




( 27 ) 
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or 


n AV 

J, = - \ + 3) ^28) 

27Tit 

where Vg is the flow velocity at the skinnner. 

Because the mass spectrometer detector used in our apparatus Is a 
density sensitive device, we are concerned with detector signals which 
are proportional to number density# In terms of number density at the 
detector, n^, equation (28) can be expressed as 

n = ^ (ym2 + 3) (29) 

2it)1 

Knuth (ref. 25) and Anderson et al. (ref. 3A) recognized that beam 
divergence would effect beam intensity at the detector and derived equa- 
tions to account for divergence. Their analysis again assumes continuum 
flow to the skimmer and collisionless flow downsfrcetm from the skimmer. 

For the usually encountered experimental case where the orifice to skimmer 
distance, Xg, is much smaller than the skimmer to detector distance, 

Knuth derives the distrlbucion 

d^ " ^ (t) ”s ® ^ T ^ cos s + (30) 

where Sg = Mg/y/i and 3 is the half angle subtended by the skimmer 
at the source as shown in figure 7. When equation (30) is integrated over 
all velocities, u = 0 to u = and over 8 from 0 to 8, a general 
equation for the number density at the detector results 


- 2Sg cos^8 j (1 + etf Sg cos 8) exp(-Sg sin^8) 


+ — (1 - cos B)S„ exp(-sj) + (2S^ - 1)(1 + erf S^) (31) 

s S S Sj 

For small skimmer openings and high Mach numbers at the skimmer, i#e,, 

2 

cos 3 i 1, S cos 3 1 and 2S 1, equation (31) can be reduced to 

s s 

the more tractable approximate form 

1 - cos ^8 exp(-s| sin^8) 




(32) 
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If parallel flow to the skinner is assuned (i.e., if divergence is 
neglected) so that sin S ~ 6 and cos B « 1, equation (30) can be inte- 
grated to give 

"d - [■^ W ^ ^ *s> * <M> 

Equations siodlar to (32) and (33) have been derived by LeRoy and 
Covers (ref. 37) , Hagena and HDrton (ref« 38) « Anderson (ref« 19) , and 
French (ref* 39) • The above equations can be ised co calculate an **ideal- 
ised** theoretical density at the detector for the case where the skiner 
is placed close enough to the source orifice to be in the continuum flow 
field. In practice hoiwever,. this case is usually avoided for high pres- 
sure stapling systems. Generally the skiner is positioned relatively 
far from the orifice to provide adequate pumping in this region and to 
reduce skimmer interference. (Usually thc' experimenter measures intensity 
as a function of to find where the experimental maximum occurs.) 

Experimental considerations thus dictate the necessity for having shiner 
to orifice distances sufficiently large to preclude continuum flow to the 
skimmer. Therefore we must now examine the case where the transition 
from continuum flow takes place upstream of the skinmer. 


(b) Transition Flow Upstream of the Skimmer 

Downstream from the onset of translational freezing a relatively 
broad velocity distribution is found in the radial direction » and this 
greatly influences the axial beam intensity at the detector. Experimental 
determinations of Tj^ show a dependence on orifice to skinmr distance 
which is intermediate between the collision dominated form 

*» Tj I « and that expected in the free zaolecular limit 

The variation of Tj^ with downstream distance from the orifice has been 
attributed to the removal of collision partners from the beam core caused 
by the relatively large transverse velocities imparted by collisions after 
the onset of translational freezing (ref. 40) > 

Hamel and Willis (ref. 41) presented a physical explanation of the 
situation which is useful to follow. As the expands decreases 

due to the geometrical effect of the source flow spherical expansion but 
in the collision dominated range energy is transferred between the par- 
allel and perpendicular modes of motion because T| | > Tj . The collisional 
drain of energy cools T|| so in the collision dominated limit T| | « 
and these temperatures vary in proportion to the reciprocal of the distance 
from the orifice to the 4/3 power. As collisions become less frequent 
still decreases due to geometric expansion but collisional cooling of T| | 
decreases and eventually freezes. A relatively small collisional transfer 
of energy still takes place after translational freezin^^ occurs but be- 
cause T|j >> Tj this has a small effect on T|i but a discernible effect 
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on the perpendicular tei^>erature and does not decay as espected in 

the free aolecular case. Thus collisions are aore effective in the rare- 
fied regine in inhibiting the decay of Therefore » does not de- 

crease as the reciprocal of distance squared as expected in the free 
aolecular limit. 

From the preceding discussion ve can see that the beam core does not 
becooe collisionless and this effects the intensity at the detector. 
Shama et al. (ref. 42) and LeRoy et al. (ref. 40) give the expression 
for the density at the detector as 



where is the perpendicular speed ratio at the skisner, i.e., the 

flow speed divided by the most probable randoa speed perpendicular to the 
flow velocity at the skinner. In the derivation of equation (34) it was 
assum ed th at sin 0 » 3 and cos 3 » Reaeabering that 

S * M/y 72 and approximating 8 ** r /Xg, where r^ is the radius of the 
skianer openings equation (34) cm be written as 



Me shall use this equation to calculate idealized theoretical densities 
because it is relatively easy to apply and because the conparisons made 
by LeRoy and Covers (ref. 37) indicate that it is a reasonable approxiaa- 
tion for our purposes. 

When using equation (35) we use two different methods of determining 
Mj^. The first method is based on a sudden freeze approximation while the 
s^ond method is based on expressions derived by Sharma et al. (ref. 42) 
froa correlations with experimental data. The first method, when applied 
to our experimental conditions, yields theoretical densities which are 
about three times greater than our eszperiaentally measured densities. 
However, the second method predicts theoretical densities which are In 
good agreement with our experimental meas -rements. 


(1) Sudden Freeze Approximation for Mi 

s 

For this method the expansion is assumed to be isentropic to some 
distance, Xq, do%mst ream from the orifice. When the orifice to skimmer 
distance, X^, is greater than Xq, the Mach number at Xq is taken as 
the terminal Mach number, i.e., M^ « Mq, and 
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M 



(36) 


This use of M]> is equivalent to assuoing continuum flow upstream of 
aud molecular flow downstream from this point. With this assumption it 
follows that • (Z-/Xg)^, idiere Tq is the translational ten^ra- 

ture at calculated from isentropic relations, and 



At X the density n 


is calculated from 


n 


q 


1 + 


(y - 1) 
2 



(37) 


(38) 


aid the density at the skiraar is assumed to decrease according to 


n 

s 



(39) 


Therefore, the idealized density at the detector can be calculated for a 
given system geometry and source gas conditions, i.e., given X^, i, Dq, 
tg, Y» Oqi ?o* ^0* 


(2) Experimental Correlations for 


M 


is 


Sharma et al. (ref. 42) correlate experimental results for monatomic 
gases to obtain equations for calculating . The data on perpendicular 

tesqueratures indicate that this temperature decreases as though the expan- 
sion were isentropic to a transition point located several times as far 
from the source as is the point at which T|; begins to deviate for an 
isentropic expansion. The distance downstream from the source orifice to 
this transition point is given by 


JLl 


D 


0 



(40) 


fdiere Tq is the characteristic time for translational relaxation of the 
gas species as given by 


^0 



(41) 
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Here a is the collision diameter and M is the molecular weight of the 
source gas at temperature Tq and with density nn. From the argon Jet 
data of Habets et al. (ref. 43), who found Sj^ « X/R (where R is a func* 
tion of source conditions), Sharma et al. develop the expression 



(42) 


For the second method calculations equation (42) is used to obtain the 
value of used in equation (35) and n^ is calculated exactly the 

same as for method one. 


APPARATUS 

Our high pressure mass spectrometric sampling apparatus is shown in 
figure 8. Before proceeding with a detailed description it is instructive 
to consider the entire system briefly in general terms. 

A schematic cross-sectional view of the beam path portion of the sam- 
pler is shown in figure 9. The atmospheric pressure gas to be sampled 
enters the system through the sampling orifice located at the bottom of 
the vertical configuration. For sampling orifice Knudsen number less 
than 10*^ the gas expands into stage I as a free jet. Stage 1 vacuum 
pumping maintains pressure ratios greater than 10^ across the sampling 
orifice. The conical skimmer intercepts the gas flow and passes the cen- 
tral portion or core of the flow into stage II where the core is periodi- 
cally interrupted by a toothed rotating disk beam chopper. Pressure in 
stage II is less than 10^^ torr. The chopped gas beam enters the stage III 
differential pumping vacuum chamber through a collimating aperture. Pres- 
sure in stage III is less than 10*^ torr The beam exits from stage III 
through another aperture and passes into stage IV and the ionizer of the 
quadrupole mass spectrometer. Pressure in stage IV is less than 10*^ torr. 
Ions formed in the electron bombardment ionizer are electrically focused 
and transmitted to the quadrupole mass filter. Selected ions transmitted 
by the filter are electrically deflected into the off-axis electron multi- 
plier whose output is electronically monitored. Thus a gas sample at high 
pressure is continuously transferred to the low pressure mass spectrom- 
eter where the gas is analyzed on the basis of the mass -to -charge ratio of 
its constituents. Maintaining the fidelity of the sample depends on spe- 
cific design parameters of the apparatus. These parameters will be con- 
sidered in detail for the various components of the apparatus. 


Vacuum Systems 

Plan views of the stages I through IV vacuum chambers are shown in 
figure 10 where the letters designate flanged ports. This portion of the 
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apparatus was designed in consultation with Drs. F. T« Greene and T. A« 
Milne of the Midwest Research Institute. The vacuum chambers » skimmer, 
beam collimator, and chopper assembly were constructed under contract by 
the Midwest Research Institute. 

All vacuum chambers are constructed of stainless steel, mechanical 
pump lines are copper, 0- rings are Vlton, ionization gage tubes are 
Bayard-Alpert type and diffusion pusq> fluid is Santovac 5. The vacuum 
system for each stage is controlled by a manual/automatic electrical con- 
trol system. In the automatic mode of operation an ionization gage is 
the primary element controlling a pneumatic high vacuum valve which iso- 
lates the vacuum chamber from the respective diffusion pump. A thermo- 
couple gage tube, located at the diffusion pump esdhaust, controls the dif- 
fusion pump power and a pneumatic valve in the fore line. Diffusion pump 
temperature and water flow interlocks are also part of the electrical con- 
trol system. 

The overall size of the vacutim system was made to make the beam path 
as short as possible while providing sufficient pumping speed for each 
stage* The stage I punq>ing requirements were most critical and thus this 
stage set the general scale of the apparatus. The size of the stage I 
vacuum chamber was dictated by conductance requirements and physical con- 
siderations related to pump size and location. Stage I is pumped by two 
liquid nitrogen trapped 10-inch diffusion puoq>s. Each trap-pump com- 
bination is attached to tlie stage I chamber, with 0-ring seals, at 
flanges G. Nominal conductance for each cold trap is 4200 liters per 
second and the maximum air pumping speed of each diffusion pump is 
5300 liters per second. Each diffusion pump is backed by a direct drive 
two-stage rotary mechanical vacuum pump with a maximum pumping speed of 
61 liters per second. These mechanical pumps are also used to rough 
pump stage I through 4-inch lines connected to the chamber by way of 
manually operated slide valves attached to flanges B. The cold trap- 
diffusion pump can be isolated from the chamber, approximately at flange 
G, by specially designed air operated lift valves. These valves' pneu- 
matic operators are attached at flanges C. Stage I pressure is measured 
with either the ionization gage tube connected to flange A, or a capaci- 
tance type gage located at flange A^« The capacitance gage has a range 
to one torr full scale and a resolution of one part in 10^. Flange E 
accepts the sampling orifice inlet assembly which will be described later. 

The stage II vacuum chamber is the rectangular parallelepiped located 
in the top-center section of the stage I chamber. The stage II chamber is 
evacuated through a 6-inch diameter pipe which physically passes through 
the stage I chamber, but is isolated from it for vacuum purposes > This 
pipe terminates at flange F. A 6-inch diffusion pun^-liquid nitrogen 
cold trap combination is connected to flange F through a 6 inch, pneu- 
matically operated, slide valve. Conductance of the trap is 1700 liters 
per second and the maximum air pumping speed of the diffusion pump is 
2400 liters per second. An 8.5 liters per second two-stage rctary mechan- 
ical pump is used to back the diffusion pump and rough pump stage 11 
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throu^ appropriate valving attached to flange A 2 « Pressure in stage II 
is loeasured with an ionization gage tube connected to flange The 

stage I and II vacuum chambers are interconnected throu^ the skimmer 
assembly which is mounted on the internal flange S shown in figure 9. 

The stage III vacuum chandler consists of a vertical cylindrical sec* 
tion and horizontal length of pipe extending to the rear and terminating 
at flange The circular cross section of the pipe is compressed into 

an elliptical section where it is welded to the vertical cylinder* 

Stages III and II are interconnected through the collimator assesd>ly 
mounted on the internal flange T shown in figure 9* A 4- inch diffusion 
pump*liquid nitrogen cold trap cond)ination is used to evacuate stage III. 
These are attached at flange I through an elbow and pneumatically operated 
slide valve. Nominal conductance of the trap is 275 liters per second and 
the maximum air pumping speed of the diffusion pump is 700 liters per sec* 
ond. An 8.5 liters per second two*stage rotary mechanical pump is used 
to back the diffusion and rough pump stage III through a pipe appropriately 
valved and connected to the elbow above the slide valve. Pressures in this 
stage are measured with an ionization gage tube connected to flange A 3 * 

The stage III vacuum chamber can be isolated from stag^ IV by a spe* 
daily designed valve. This bellows sealed , manually operated valve is 
located in stage III. Valve opening and closing are accomplished by 
translating a ball mechanism. In the open position the ball is retracted 
to provide a clear passage from stage III to stage IV. In the closed posi- 
tion the ball seats on an 0*ring installed around a 0.32 centimeter diam* 
eter opening in flange H. The valve mechanism is operated through a bel* 
lows attached to flange V. 

The cylindrical tee*section attached to flange H comprises the stage 
IV vacuum chamber. This chamber houses the quadrupole mass spectrometer. 
All vacuum connections to stage IV are made with copper gasket seals. A 
1000 liters per second (air) baffled getter*ion pump, attached to flange J 
through an elbow, is used to evacuate stage IV. Hie mass spectrometer is 
mounted on a flange which mates with flange K. Pressure is measured with 
an ionization gage tube connected to flange Ag. 


Sampling Orifice Inlet Assembly 

The three assemblies shown in figure 11 have been used for various 
experiments. Flange L, which mates directly with flange £ of the stage I 
vacuum chamber, and spacer ring LS is common to each assembly with either 
flange M 2 , or M^. Spacer rings of various height are used to give the 
desired orifice to skimmer distance. Four screws pass radially, 90 degrees 
apart, through flange L as shown in figure 11. These screws facilitate 
lateral movement of the orifice flange, M, for alinement purposes. 

Flanges or M 3 are used when it is desired to allow the sampling 
orifice to get hot, and flange M 2 is used when a cooled orifice is desired. 
In the Mj configuration the platinum orifice cone is welded directly to 
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the stainless steel plate • Here n type I orifice is used and its geometry 
is given in figure 12. fox the M 2 configuration a type II orifice is 
soldered directly to the water cooled plate. The M 3 configuration uses a 
type II orifice welded to the end of the stainless steel tube ertension 
attached to the stainless steel plate. In configurations and M 3 the 
stainless steel plate is attached to a water cooled plate and vacuum seal 
is made with an 0«-ring. The water cooling in these configurations is jsed 
only to protect the 0 -rings. 

With the two types of platinum-rhodium alloy sampling cones shown in 
figure 12 9 the orifice is located in the ^ex of the cone. The orifice 
is machined so as to produce a sharp edged circular opening with zero 
channel length. 


Vertically Translatable Orifice Inlet Assembly 

For experiments where it was necessary to have a continuously vari- 
able orifice*t:o-skimcier spacing and reduced source gas pressures » the spe- 
cial assembly shown schematically in figure 13 was used. This unit was 
added to the apparatus by substituting it for flange M of figure 11. By 
turning the screw collar the orifice to skimmer distance could be continu- 
ously varied with a precision of better than 0.02 centimeter. The cylin- 
drical source cdiaai>er is keyed to the housing so that the chamber and ori- 
fice do not rotate during vertical translation. 

Pressure In the source chand>er could be varied from a few microns to 
800 torr by suitable vacuum pumping and pressurization from gas cylinders. 
A mechanical pump was used for evacuation of this chamber. Pressures be- 
low 1 torr were measured with a capacitance type vacuum gage. For source 
chamber pressures above 1 torr a precision dial type gage was used to 
measure pressure. 

When sampling from the source chanber, the gas under investigation 
was continuously admitted to the chamber through leak valves and the flow 
rate was adjusted to maintain the desired pressure in the source chamber. 

A type II sampling cone was soldered to the upper end of the source cham- 
ber and the orifice diameter was 0.025 centimeter. The general design of 
this unit provided good stage I pumping in the vicinity of the orifice 
even at small orifice^to-skimmer distances. 


Skimmer and Ck>lliinator Assemblies 

Cross sectional schematic views of the skimmer and collimator are 
shown in figures 14 and 15, respectively. The flange mounting portion 
of these assemblies are similar and they will be described together. The 
skimmer and collimator flanges mate with an intermediate flange as shown 
in figure 14. The intermediate flanges in turn mate with flange S (fig. 9) 
for the skimmer and flange T (fig. 9) for the collimaLor. The intermediate 
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flange in each case is located with dowel pins to the flanges S or T« A 
screw retaining ring serves as a clamp. Vhen the retaining screws are 
secured in place and the lockdown screws are removed, the assembly can be 
removed from the apparatus. The assembly can then be replaced in the 
apparatus, using the dowel locating pins, in exactly the same place as 
before removal. With both the retaining and lockdown screws in place but 
not tightened, the skimmer or collimator can be moved laterally for aline- 
ment purposes. Once alignment is achieved, securing the retaining screws 
fixes the position of the skimmer or collimator, respectively. 

The skimmer cone is span from 0.025 inch thick copper. The opening 
in the apex of the 60 degree cone is machined in a way to provide a sharp 
edged circular opening with zero channel length. For all the experiments 
reported here a single skimmer was employed and it had a 0.076 centimeter 
diameter opening. 

The collimator cone is made similarly to the skimmer but with an 
included angle of 120 degrees and a 0.159 centimeter diameter opening. 


Beam Chopper Assembly 

The beam chopper assembly consists of a motor driven segmented disc, 
a light bulb and a photo diode detector. This assembly is located in 
the stage II vacuum chamber and attached to the collimator as shown sche- 
matically in figure 16. Th«. rotating disc is made of 0.032 inch thick 
aluminum. The disc is shaped as shown in figure 17 and carefully balanced 
to yield vibration free rotation. 

A chopping frequency of 150 hertz was routinely employed in our 
studies. The motor rotation speed was electrically adjusted to produce 
the desired frequency as measured with the lamp-diode combination. Elec- 
trical connections to the beam chopper were made through an insulated 
feed through attached to flange A7 shown in figure 10. 


Quadrupole Mass Spectrometer 

An Extranuclear Laboratories Inc. quadrupole mass spectrometer system 
is used in our apparatus. This system is composed of the following model 
number components: 041-1 axial ionizer, 021-1 electron energy-emission 

control, 4-270-9 ELFS quadrupole mass filter with 1.6 ^ 22 0 cm poles, 

011-15 power supply - 1 to 6 MHz, 1.25 MHz high-Q head, 031-3 electrometer, 
032-5 preamplifier, channeltron electron multiplier and 091-6 digital mass 
programmer. The quadrupole mass filter has a nominal range from 0 to 
500 AMU but in practice we have been able to operate above 600 AMU. 

The entire spectrometer is mounted on an 8-inch diameter copper gasket 
type flange which mates directly to flange K of the stage IV vacuum cham- 
ber. All electrical connections to the spectrometer are brought out through 
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its mounting flange. This flange Is also equipped with a window to pro- 
vide a clear line of si|^t through the entire spectrometer. An electri- 
cally operated shutter protects the window from condensible species in 
the beam. 

Entrance into the mass spectrometer electron bombardment ionizer is 
through a 0.32 centimeter diameter aperture. The mass spectrometer is 
physically positioned in the stage IV vacuum chamber to locate the en- 
trance aperture approximately 1 centimeter above the stage III exit aper- 
ture. Three bellows sealed micrometer screws are located at the lower 
end of the stage IV chamber 120 degrees apart radially. These screws 
bear against the quadrupole mass filter pole housing and facilitate cen- 
tering the entrance aperture with respect to the stage III exit opening. 
This positioning can be accomplished while stage IV is under high vacuum. 


Signal Processing 

Beam modulatlcm and phase sensitive detection techniques are used in 
signal processing to Increase the signal -to-noise ratio, S/N. The signal 
of Interest is the ion current produced by ionization of the chopped sam- 
ple gas beam. There are numerous sources of noise signal and one of 
these is the signal due to background gas which finds its way into the 
ionizer. 

Beam modulation is achieved with the chopper assembly which both 
periodically interrupts the gas beam and produces a square wave refer- 
ence voltage at the modulation frequency. A schematic wiring diagram of 
the chopper system is shown in figure 18. The variable frequency oscil- 
lator output is adjusted in frequency and amplitude to make the motor 
speed of rotation appropriate to result in the desired chopping frequency. 
A regulated fixed current is passed through the lanq> and as the sectored 
disk allows light to illuminate the photo diode, an output reference sig- 
nal voltage is produced. This square wave reference signal at the modu- 
lation frequency is used to drive the synchronous gate in the phase sensi- 
tive detector signal processor. A Princeton Applied Research model HR-8 
lock-in amplifier is the signal processor. 

The output of the spectrometer electron multiplier is fed directly 
into the 032-5 preamplifier mounted in close proximity to the signal feed 
through in the spectrometer support flange. The preamplifier contains 
remotely selectable input resistors of 10^, 10^, and 10^ ohms. The out- 
put of the preamplifier is routed in parallel to the electrometer, the 
lock-in amplifier, an oscilloscope, and the digital mass programmer. Sig- 
nal routing is shown in figure 19. 

The d.c. output of the electrometer is proportional to the total ion 
current, l.e., the sum of (1) the average value of the ion current pro- 
duced from the chopped gas beam, (2) the ion current resulting from back- 
ground gas which enters the ionizer, and (3) noise current. The elect^om- 
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eter output can be filtered with preselectable time constants to smooth 
the signal. 

The d.c. output of the lock- In auq)llfler is proportional to the compo- 
nent of the total Ion current that Is In phase with the modulated gas beam 
pulses. Of course the lock-in amplifier must be appropriately tuned with 
regard to synchronization frequency and phase angle. The lock-ln amplifier 
output can also be filtered with preselectable time constants. 

Electrometer and lock-in amplifier outputs are recorded on separate 
channels of a two channel high speed servo recorder. 

The quadrupole mass filter is controlled by a power supply which pro- 
vides three options for making the filter sweep through an AMU range. In 
the manual mode, the operator manually adjusts the control potentiometer 
to sweep the mass filter. A fixed setting of this potentiometer corres- 
ponds to a fixed AMU position for the filter and thus one can tune the 
filter to rest on any AMU of Interest. In the sweep mode the filter sweep 
rate follows the saw tooth ramp rate selected for the oscilloscope. *^n 
this mode the AMU starting point and range are Independently selectable 
by the operator. The third sweep mode Is labeled external command. In 
this mode the sweep rate, starting AMU, and AMU range are controlled by 
commands received from the 091-6 mass programmer. The operator digitally 
programs the mass programmer; there are sixteen channels which can be 
programmed. A single channel can be made controlling or all sixteen channels 
can be run through sequentially. 

The mass programmer also contains a signal processing section where 
the ion current signal can be amplified by a factor from 0.01 to 99.9 
according to programmable command for each of the sixteen channels. This 
has proven to be a useful feature which we alternately employ by suitable 
switching as shown in figure 19. 


PROCEDURE 

Alignment 

Optimum performance of the sampling system is critically dependent 
upon proper alignment of all orifices and apertures in the gas beam path. 
A low power (0.2 mW) helium-neon laser is used to align our system. The 
laser is supported by a mechanism which facilitates setting the beam ver- 
tically and translating the beam in two directions perpendicular to the 
vertical. The exit opening from stage III is fixed in position on the 
apparatus and therefore all alignment is done relative to this element. 
The laser is pointed up through the system from the bottom and the beam 
is centered in the stage III exit opening with all other elements removed 
from the apparatus. Next the quadrupole mass spectrometer is secured in 
place, and the three centering micrometer screws in the stage IV vacuum 
chamber are adjusted until the laser beam is seen to emerge clearly 
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through the window at the upper end of the spectrometer. At this point 
the ionizer inlet aperture and stage III exit aperture are lined up cor- 
rectly* The collimator assembly is then positioned and adjusted to pass 
the laser beam. Similarly the skimmer asseably is added and brought into 
alignment* Finally^ the sampling orifice assembly is put in place and 
positioned to pass the laser beam* 

When the sampling apparatus is pumped down and put into operation, 
the sampling orifice position is given a final ''fine tuning" adjustment. 
This adjustment is made while ambient laboratory air is being sampled. 

The mass filter is set on the maximum of the nitrogen peak and the elec- 
trometer output on the recorder is observed while adjusting the sampling 
orifice centering screws in flange L. Adjustments are made to maximize 
the nitrogen signal. 


Spectrometer Tuning 

The quadrupole mass spectrometer used in our apparatus is one for 
which all settings of the ionizing, focusing, and mass filtering elements 
must be determined by the experimenter. The detection of low level sig- 
nals depends very strongly upon the adjustment of these elements. Cptl- 
mum adjustment corresponds to maximizing the signal- to- noise ratio (S/N) 
rather than just maximizing the signal. 

Our experience has indicated that the two controls which most strongly 
affect maximization of S/N are the emission current and the ion energy* 
Although in general the signal strength increases as each of these param- 
eters is increased, a point is reached where the noise level begins to in- 
crease more rapidly than the signal level. We have found that optimum 
S/N is achieved at relatively low emission current (typically about 3 mA) 
and ion energy (typically 10 V at 32 AMU) . 

Tuning is accomplished by observing an intense peak on the oscillo- 
scope and adjusting the controls to produce the proper peak shape. This 
is usually done by observing the nitrogen peak while saiopling ambient 
laboratory air. A fast repetitive sweep rate is used and the resolution 
control is set to a low value so that the peak covers several AMU and is 
about 8 centimeters wide on the oscilloscope screen. The tuning controls 
are adjusted so that the top of the peak is flat, the high mass side of 
the peak is vertical, and the low mass side is a straight line with maxi- 
mum slope* These features of the peak shape are sought while also endeav- 
oring to maintain as large a peak amplitude as possible. The ideal shape 
is shown in figure 20 by the heavy solid line. Once the tuning is opti- 
mized, the resolution can be increased considerably before the peak ampli- 
tude is affected* As the resolution is increased, the position of the low 
mass side of the peak should move toward higher mass as shown by the dashed 
lines in figure 20. When the low mass side ol the peak just intercepts the 
high mass side, e.g*, point A in figure 20, the spectrometer is tuned to 
the highest resolution which gives the maximum S/N. Increasing the reso- 
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lution control f«tthtr Increases the resolution but only at the expense 
of a loss In sensitivity, e.g. , point B in figure 20. 

Optimum resolution and sensitivity at low AHU does not necessarily 
represent the optimum at high AMU. The control which most strongly af- 
fects the ratio of the low to high AMU resolution is the AM control. 

We have found that for our spectrometer a setting of AM » 0 should be 
used to obtain resolutions which are equally satisfactory at low and high 
AMU while maintaining good sensitivity. 

The intensity and resolution are usually checked daily by observing a 
one atmosphere oxygen beam. The oxygen used contains a low level xenon 
impurity which provides a suitable signal at relatively high AMU. The 
oxygen peak and one of the xenon peaks are alternately observed while fine 
tuning. 

Typical resolutions (defined using the full-width at half -maximum) » 
used when detecting low level signals, are IX to 2X the mass number in 
the AMU region near xenon and 2X the mass number in the region near oxygen. 


Beam Modulation 

Fite (ref. 44) has pointed out that in selecting a chopping frequency 
for performing modulated beam mass spect rometric measurements two physical 
considerations govern. First, the frequency should be made sufficiently 
high to insure that the Fourier component of the background gas pressure 
fluctuations is negligibly small. Second, the frequency should be made 
sufficiently low to insure that coherence of the modulated beam is not 
destroyed as the result of velocity spreading of the beam pulses. 

For effusive molecular beams where there is a broad distribution of 
velocities, the second consideration above is of critical impottance. 

Fite and Brackman (ref. 45), Harrison et al . (refs. 46 and 47), and 
Yomamoto and Stickney (ref^ 48) have analyzed the effusive beam case and 
worked out the dependence of the amplitude and phase shift on frequency, 
molecular mass and flight path length. 

For supersonic beams. Miller (ref. 49) has shown that cue second con- 
sideration above is not critical. Miller showed that amplitude correc- 
tions are minor when the Mach number i*^ greater than about 5. This 
results from the fact that at high Mach numbers the velocity distribution 
is sufficiently narrow that rhe spread in velocities does not effect the 
coherence of the beam even at relatively high chopping frequencies, e g , 
1400 hertz. 

Because our sampling system produces a high Mach number supersonic 
beam, we have only considered background gas pressure fluctuation effects 
in selecting a chopping frequency. We have used the analysis ot Fite 
(ref. 44) to estimate the ratio of amplitudes of pressure fluctuations, 
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in the stage II vacuum chamber, for a chopped beam as compared to an un- 
chopped beam. The stage II vacuum system's tijpe constant is approximately 
40 milliseconds. At a chopping frequency of 150 hertz the pressure fluc- 
tuation amplitude is attenuated by a factor ot about 40 and at a fre- 
quency of 400 hertz the attenuation factor is 100. 


RESULTS AND DISCUSSION 

The general performance of our sampling system was irltially assessed 
by sampling the ambient laboratory atmosphere throufh a 0,025 centimeter 
diameter orifice. The orifice-to-skimmer distance was arbitrarily set at 
Xq/Dq “ 50. These initial tests demonstrated that the vacuum pumping sys- 
tems met the design criteria and were adequate for ambient temperature, 
atmospheric pressure sampling through a 0.025 centimeter diameter orifice. 
The pressure in stage I was about 1.5 microns, and this value corresponds 
to the region where the lO-inch diffusion pumps begin to achieve their 
maximum pumping speed. This result indicates that the 0,025 centimeter 
diameter orifice is about the maximum that can be used on our sampler with 

the stage I pumping provided. The pressure in stages II, III, and IV was 

2x 10"5^ 3>«.10*^, and ^10*® torr, respectively. 

A portion of a typically recorded spectrum for the laborator> ambient 

atmosphere is shown in figure 21. The peaks at 32, 33, 34 AMU are due to 

oxygen and those at 36, 38, and 40 AMU are due to argon which exists in 
the atmosphere at about 0.93 percent by volume. An estimate of the sam- 
pling system's sensitivity was ascertained by observing that the Ar^^ iso- 
tope was readily identified and measured here with the spectrometer and 
electronics set to only an Intermediate sensitivity levels This isotope 

is of the order of 5 parts per million (ppm) in the ambient atmosphere. 

The 37 AMU peak in figure 21 was attributed to the ion species 

H(H20)2* a subsequent scan of ambient air was made to higher AMU and is 

shown in figure 22. This scan was at a higher sweep rate and lower sensi- 
tivity than that shown in figure 21, The water and water cluster peaks 

are clearly visible and in accord with observations made by other investi- 
gators (ref. 50), To confirm these observations we saoipled ambient air 
which was first passed through a drying column. A scan ot the dry air 

sample is shown in figure 23, The water peak at 18 AMU is eeen to be 

greatly reduced and the water cluster peaks are absent. The scans shown 
in figures 22 and 23 were made at the same sensitivity, scan speed and 
time constant. At higher sensitivity, slower scan speed and increased 
time constant the dried air sample still showed some water clustering 
peaks and this demonstrated that, as expected, the drying v; as not complete. 

After the preliminary testing described above, experiments were per- 
formed to determine (1) how measured intensity varied as a function of 
Xq/Dq and source pressure, Pq, (2) how the measured llitenslty of an argon 
beam compared to theorecically calculated intensity, and (3) the ultimate 
sensitivity of the sampling system. For these experiments the transljf- 
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able sampling orifice assembly was Installed on the sampler Experiments 
were performed with argon, nitrogen, and oxygen as the room temperature 
source gas. 


Argon 

Argon was sampled at various source pressures and the iii*phase ion 
current was measured as a function X^/Dq. Typical results at three 
source pressures are presented in figure 24. The Intensity of the argon 
signal is seen to pass through a rather broad maximum near X^/Dq » 120 
and then fall sharply at lower values of X^/Dq. The sharp drop in in- 
tensity at lower orifice to skimmer distances Is similar to that reported 
by other Investigators (ref. 51>, and this drop is usually attributed to 
skimmer interference and shock effects (refs. 51 and 52). We believe 
that in our case the drop in intensity at low values of Xg/Dg may also 
be due to background gas scattering effects because at small orlfice^to^ 
skimmer distances the pumping in the region of the orifice is probably 
somewhat restricted. The gradual decrease in intensity at Xg/DQ values 
beyond the maximum may be due to the decreased density at the skimmer 
entrance resulting from ti.e increasing ori f ice-to-sklmmer separation 
According to the idealized theoretical analysis of Parker et al. \xef. 36), 
the flux at the detector varies as (Xg/DQ)~2/3 for argon. This is in 
reasonable agreement with our observations for Xg/Dp values greater 
than where the maximum intensity occurs 

It is Interesting to note that if one makes a simplistic calculation 
of a skimmer Knudsen number at the point of maximum intensity, a value 
is obtained which is about equal to that usually considered to be a cri- 
teria for free molecular flow. For example, consider the case cf argon 
with Pq » 700 torr and Xq/Dq = 120 As will be shown later, the num- 
ber density at the skimmer is calculated to be = 2 35 00^^ aicms-cra“3. 
The mean free path at the skimn^er, is ^aiculeted to be U- 7i centimeter 

from the expression where fhe collision cross section 

0 is taken to have the value 3.67-10*^ cm* For the 0,076 centimeter 
skimmer diameter, the skimmer Knudsen number, Kn^ , is calcul.ted from 
9,3, Because a value of Kng = 10 is taken as the 
criteria tor free molecular flow, it is reasonable to conclude that for 
our example free molecular flow prevailed where maximum intensity was 
achieved. Upstream of where free molecular flow prevails, collisions 
still play a role and shock or scattering effects would be expected to 
attenuate the beam as previously postulated. Thus this qualitative pic- 
ture in general can explain fhe experimentally observed shapes of the in- 
tensity versus Xg/D|Q curves. Furthermore, this approach is in general 
accord with the findings of other investigators who observed similar curve 
shapes (refs. 53 and 54). For example, Scott and Drewry (ref. 54) con- 
cluded that experimental results rise to a maximum and tend to approach 
theoretical predictions as Kn^ increases. 
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Our aeasured dependence of intensity on source pressure for 
Xg/I^ s 120 is shoim in figure 25. For pressures below about 600 tort 
the intensity is seen to be directly proportional to source pressure and 
this is in accord with theoretical considerations. The fall-cff of in* 
tensity at higher pressures was not expected. We speculated that this 
inigjht be related to stage 1 pumping. At these higher pressures the 
stage I pressure is in the region where the diffusion pumps just begin to 
attain their maximum pumping speed. However , the stage 1 pressure is 
almost linearly related to source pressure, as shown in figure 26, even 
to source pressures of 800 torr. One would expect to see a departure 
from linearity at higher pressures if the diffusion pumps were not pump- 
ing at their maximum speed. 

Pressures in stages II, III, and IV were measured as a function of 
source pressure and conditions the stage IV pressure 

was always less than 10*^ torr and this fact demonstrates that even for 
argon the stage IV pumping capacity was more than adequate. The pressures 
measured in stages II and III were found to be a function of source pres* 
sure and ^ shown in figure 21 where results are presented for 

two source pressures. 

To compare measured intensities with theoretically calculated values 
it is necessary to relate measured ion current to beam flux or number 
density which can be expressed in terms of source conditions. Because 
the ion source of the mass spectrometer is a density sensitive detector, 
we chose to characterize the beam in terms of number density n. The 
cent^itline number density at the ion source, n^, is proportional to the 
measured ion current, I, 


* Cl (A3) 

where C Is the proportionality constant and the current is the in-phase 
component of the modulated beam current. The proportionality constant 
may be a function of molecular species but for any given species, such 
as argon, the proportionality should be equally applicable to an effusive 
or supersonic beam. If the effusive beam is subscripted by £ and the 
supersonic (nozzle) beam is subscripted by N, we can write for a partic* 
ular species i , 


.C - 

I 




i^E 





(44) 


Effusive argon beams were produced by sampling argon source gas at 
pressures in the micron range and the results are shown in figure 28 
(^nerally, effusive flow is considered to exist for orifice Knudsen num- 
bers above about 10 but out data was found to be very linear to source 
pressures near 400 microns where the Knudsen number is approximately 0.5. 
The data to Pq = 400 microns was fitted to a ttc *ght line by the least 
squares method and the line shown through the data )oints in figure 28 
was found to fit the equation. 
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I_ » l.95xlO"^^ . p- - 1.04«10”^ (45) 

iiT b V 

At a source pressure of 30 nicrcns the mraber density at the detec- 
tor (ion source) was calculated to be Ar^d^ * 1*81 »10^ oolecules ca*^ 

for 1^/Dq * 120 and a skisnier^o^detector distance, t, equal to 
42.82 centiaeters. Fron the least squares line equation a source pres- 
sure of 30 aicrons results in a beam current of 3.75>^10^^^ aiq>ere« The 
proportionality constant At^ i^ then calculated to be Ar^ “ 3.15 ^-lO^^ 
molecules cmT^ ampere*^. 

Supersonic argon beams resulting from a source pressure* of 700 torr 
were measured as a function of X3/DQ and the measured ion currents, 

Afig, were converted to number density, Ar^d^» multiplying Ar^N 
ArC* Results are presented in figure 29 where it can be seen tlmt the 
maxi mum density occurs at X^/Dq * 120 where Ar°d^ 3.07*10^^ molecules 
cm*^. 


Theoretical densities at the detector can be calculated by equa- 
tion (35) and the techniques described in the section on Beam Formation. 
Calculations were made for the source conditions Pq - 700 torr of argon, 
Tq a 298 K, Dq = 0.025 centimeter, X3/DQ = 120, Dg * 0.038 centimeter, 
and I “ 42.82 centimeters. When Mx was calculated according to equa- 
tion (37) a ^(Xg/Xq)] and the first method, we obtained 

A^d^ ® 1.03x10^2 molecules cm"^ for a theoretical density. This value 

is a factor of about 3.4 higher than our measured value. When was 

calculated according to the second method and equation (42) 

(Mi^ = (2.50//7HXg/DoHDo'(’^0^0^~®‘^J^» obtain * 3.24-10ll 

oolecules ca~^ which is only about 5 percent higher than our oteasured 
value. It must be recognized that while both theoties are based on 
many idealized assumptions, the second method calculation results from a 
theory derived from correlations with experiment. Therefore, it is not 
surprising that the density derived from the second calculation is in 
better agreement with our measurements. 

For either method of calculation of the theoretical density at the 
detector, the agreement with experiment is considered to be good and an 
indication that at X^/Dq * 120 skinmier interference is not a serious 
problem in our apparatus. The calculations also yield additional informa- 
tion which gives some insights into the sanq>ler*s operation^ For the 
source conditions noted Mj « 30.2 and the transition to "collisionless 
flow” occurs near X^/Dq * 28 where = 4 26 10^^ molecules cm“^. The 
density at skimmer is Hq * 2. 35 >"10^^ molecules cm*^. 

A final series of argon beam experiments were performed to ascertain 
if abnormal condensation effects existed in our sampling system beam for- 
mation process. Condensation phenomenon are an important consideration 



28 


in high pressure sampling and extensive research (refs. 55 to 58 and ref- 
erences cited therein) has been devoted to this subject. VHiile most as- 
pects of condensation in free jets is beyond the scope of our discussion^ 
Milne and Greene (ref. 55) have shovn hov measurement of the argon dimer 
to monomer ratio can be used to evaluate the performance of a sampler. 

Argon source gas was sampled at room ten 9 >erature through a 0.025 cen- 
timeter diameter orifice. Orifice-to-skimmer distance was held constant 
at the optimum determined previously, i.e., X^/Dq « 120. In-phase ion 
currents were measured for the argon monomer and dimer. Measured ion cur- 
rents were not corrected for mass spectrometer mass discrimination, multi- 
plier gain or sampler mass discrimination because such corrections were 
considered to be relatively small for the short range from 40 to 80 AMU. 

At a source pressure of one atmosphere we measured an argon dimer to 
monomer ratio of 3.6x10^^. This value is in good agreement with Milne 
and Greene's data (ref. 55, fig. 5) extrapolated to our orifice size. We 
therefore conclude that for the orifice size, source pressure and tempera- 
ture, and orifice-to-skimmer distance used we had a relatively undisturbed 
free-jet expansion and beam formation. 


Nitrogen 

Nitrogen was used as a source gas to perform experiments similar to 
those done with argon. The results obtained with nitrogen were generally 
the same as those obtained with argon. Pressures in the various vacuum 
chambers as a function of source pressure and Xg/DQ were in substantial 
agreement with those measured with argon as the source gas. Intensity 
also varied with Pq and Xq/Dq in a fashion similar to argon. Effu- 
sion experiments yielded a value for the proportionality constant 
N 2 C ** 2.97x10^^ molecules cm“^ ampere”^. The measured density at the 

detector, N 2 ^d>j* plotted as a function of X^/Dq in figure 30 for a 

source pressure of 700 torr. Again the density is seen to pass through a 
broad maximum which peaks at Xq/Dq ♦ 120. The maximum density for nitro- 
gen is 1.26x10^^ molecules cm~^. 

The maximum density measured at the detector for argon beams is about 
a factor of 2.5 higher than that measured for nitrogen beams under the 
same conditions^ If one calculates a theoretical density by the methods 
of equation (35), a higher density is obtained for nitrogen compared with 
argon. The value of such calculations for nitrogen are, however, some- 
what suspect because the theoretical equations were developed mainly on 
the basis of a monatomic gas. 


Oxygen 

Preliminary experiments with oxygen as the source gas indicated that 
a xenon impurity contained in the oxygen might be useful for determining 
the ultimate sr isitlvity of our sampling system. Therefore, a number of 
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experimeats were perfonned with oxygen. Missile grade tank oxygen was 
used at atiBospheric pressure with a 0.025 centimeter sampling orifice 
diameter and with » 120. The oxygen was independently analyzed 

by conventional mass spectrometrlc techniques which had a sensitivity 
of about 1 ppm. No xenon was detected in that analysis and we therefore 
took the xenon level to be 1 ppm. 

When the oxygen was sampled with our sampling system xenon was 
readily identified. A typically recorded spectrum of the xenon isotopes 
is given in figure 31. The total current (upper trace) was obtained at 
maximum spectrometer sensitivity for acceptable resolution ^ i.e., maximum 
electron multiplier gain and maximum electrometer sensitivity. The major 
xenon isotopes are just discernible and for the 132 isotope the signal-to- 
noise ratio is about 2. h'^wever^ the in*-phase current (lower trace) 
clearly shows the xenon Isotopes and the signal-to-^noise ratio for the 
132 isotope is about 10. The in-phase current was measured with the 
lock-in amplifier sensitivity set a factor of 10 below maximum and a 
10 second time constant was used. Figure 31 demonstrates the signal^to^ 
noise ratio enhancement that is achieved by beam modulation and phase 
sensitive detection. Had a longer time constant and slower scan late 
been used, an even greater enhancement could have been attained. 

The 132 xenon isotope ±s about one fourth of the total xenon level 
and thus for the presumed xenon impurity level of 1 ppm the 132 signal 
represents about 0.25 ppm. Clearly, if the 132 signal were one fourth 
the level shown in figure 31, it would still be easily measurable at 
the sensitivities employed. On this basis we therefore expect an ulti- 
mate sensitivity for our sampling system in this general mass range to 
be on the order of 0.05 ppm. 


Other Source Gas Systems 


We have used our sampling apparatus to study reaction product forma- 
tion in flames and flowing gas systems (ref. 13) as well as to study 
volatile products emanating from metal samples subjected to various gas- 
eous environments at elevated temperatures (refs. 14 and 15). For these 
experiments atmospheric pressure gases were sanq>led through a 0.022 centi- 
meter diameter orifice in the apex of a 0.6 centimeter high sampling cone 
welded to a 2 centimeter diameter by 1.5 centimeter high stainless steel 
tube attached to the flange plate. The skimmer diameter was 0.081 centi- 
meter and 
were 1.5>10 

pheric pressure source gas. When the gas temperature was 1000*^ C, these 
pressures were 7>10*^ and IxlO*^ respectively. Stage III and IV 

pressures were less than 10“^ and 10*® torr, respectively, for both source 
gas temperatures. 


/Dq « 144. For this arrangement the stage 1 and II pressures 
and 8x10*® torr, respectively, with room temperature atmos- 


Flat flames were produced with a specially designed burner (ref. 13) 
and fuel-lean methane-oxygen flams doped with SO 2 , H 2 O, and NaCl were 
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used to study the formation of Na 2 S 04 (g)« Gaseous flame reaction 
products were measured by recording the in-phase ion current at appro- 
priate values of AMU. Composition profiles were determined by measur- 
ing each species ion current as a function of vertical distance between 
the sampling orifice and burner surface » denoted by Z. Clogging of the 
sampler orifice by condensible species was a problem which limited the 
duration of each safl^>ling experiment to about 30 minutes. 

Measured composition profiles for a 9.5 mole ratio O 2 /C &4 flame 
doped with SO 2 , H 2 O 9 and NaCl are presented in figure 32. The profiles 
observed for 02 » H 2 O, and CO 2 are similar to those reported by other 
investigators for CII 4 -O 2 flames (ref. 59). The direct observation of 
Na 2 S 04 (g), NaS 02 (g)» and NaS 03 (g) are new findings and the shape of 
their profiles is discussed in reference 13. 

The formation of Na 2 S 04 (g) was also studied in a flowing gas reac- 
tion tube system at 1413 K. The reaction tube was positioned at the 
sampler inlet in such a way that the sampling orifice on the stainless 
steel tube extended into the reaction tube and the hot zone. The orifice 
was allowed to get hot before the water-saturated oxygen, SO 2 and NaCl(g) 
were flowed through the hot reaction zone. 

Intensities were measured for the gaseous species O 2 , H 2 O, SO 2 , SO 3 , 
NaCl, Na 2 S 04 , and NaS 03 with known flows of O 2 , SO 2 , H 2 O, and known con- 
centrations of NaCl(g). The species NaOH(g) and HCl(g) with parent ions 
at m/e 40 and 36, 38 could not be measured because large in-phase back- 
ground peaks were present at these m/e values. These background peaks 
were due to argon which is a 2300 ppm impurity in the missile grade oxy- 
gen that was used. 

Equilibrium thermodynamic calculations of the coi^osition of the 
reaction products were made as a function of teiiq>erature. In the calcu- 
litions the composition of the reactants was held constant at the weight 
fractions 0.96 O 2 , 1.2x10**^ NaCl, 0.012 H 2 O, and 0.032 S 02 « Calculated 
reaction product ion intensity fractions are coiiq>ared with the equilibrium 
calculated mole fractions in the histogram shown in figure 33. 

Considering that the measured intensities were not corrected for 
sampling factors, the quality of agreement between experiment and calcu- 
lation is considered quite good. 

Our sampler has proven to be valuable in studying the atmospheric 
pressure oxidation of various metals subjected to various contaminant 
gases added to the oxygen environment (refs. 14 and 15). The high pres- 
sure mass spectrometric sampling technique has been used to identify vol- 
atile products emanating from the metal samples. For these experiments 
the metal specimen is suspended in the reaction tube and the sampling 
orifice is located 1 to 2 centimeters above the specimen. An example of 
the type of results obtained is shown in figure 34 for the system chromium 
specimen plus O 2 + NaCl(g) ± H 20 (g) at 1020^ C. With this system the new 
species (NaCD^^CrO^ with x » 1 , 2 , or 3 , (NaOH)yCrOj with v ■ i or 2 , 
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NaClHa 2 Cr 207 » and HaCBHa 2 Cr 207 were readily identified. At first glance 
one ai^t suspect that these coiqtlexes are the product of gas phase 
reactions occurring in the sai^ling process, but independent evidence 
(ref. 14) suggests that this is not the case. Even in the case of tiie 
NaCl polyners condensation clustering does not sees reasonable when one 
considers the low NaCl(g) pressure and the fact that the gas is at ele- 
vated tee^rature (ref. 60). 
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TABLE I. - MACH NUMBER EXPANSION 


EQUATION CONSTANTS 


Y 


A 

5/3 

0.075 

3.26 

7/5 

.04 

3.65 

9/7 

.85 

3.96 


TABLE II. - POLYNOMIALS FOR EXTRAPOLATING MACH NUMBER 
TO SUBSONIC FLOW REGIME 


mm 

Range 

Polynomial 

5/3 

0 <. (X/Dq) i 1.665 

M = 1.000 + 1.671(X/Dq) + 0.047.j(X/Dq)^ 

5/3 

-0.897^ i (X/Dq) i 0 

M = 1.000 + 1.671(X/Dq) - 0. 6912 (X/Dq) ^ 

7/5 

0 < (X/Dq) :< 2.228 

M = 1.000 + 1.441 (X/Dq) - 0.019q(X/Dq)^ 

7/5 

-1.041 < (X/Dq) :< 0 

M = 1.000 + 1.441 (X/Dq) - 0. 4433 (X/Dq) ^ 

9/7 

0 i (X/Dq) 2.852 

M = 1.000 + 1.174(X/Dq) - 0.015q(X/Dq)^ 

9/7 

-1.278 ^ (X/Dq) 0 

M- 1.000 + 1.174(X/Dq) - 0.239q(X/Dq)^ 
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Figure 4. - Variation o' *tressure, density, and temperature 
with dimensionless flew time for three values of the spe* 
cific best ratio y. 
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Figure 9. - Schematic cross section of high pressure mass spectrometric sampler. 






Figure' 13. - Schematic cross-sectional view of translatable orifice inlcft assembly. 
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Figure 15. - Collimator. 



Figure 16. - Beam chopper assembly. 



Figure 17. - Beam chopper segmented disc. 
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Figure 18. - Schematic wiring diagram for chopper system. 
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Figure 21. - Pcrtion ol mass spectral scan tor ambient laboratory air. 



Figure 22. - Mass spectral scan of ambient laboratory air. 
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Figure 23. • Mass spectra' scan of ambient laboratory air passed through drying column. 
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Figure 24. - Argon beam intens^y variation as a function of onfice-to-skimmer 
distance and source pressure. 
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Figure 21 . * Variation ot second and third vacuum chamber pressures with orifice to skimmer distance at two 
source pressures. 
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Figure # - Argon beam intensity as a function of source pressures m ine 
effusiw»anqe. 
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Figure 29. * Argon beam density at :.ie detector plotted versus orifice to 
skimmer distance. 
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Figure 30. - Nitrogen beam density at the detector plotted versus orifice to* 
skimmer distance. 
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Figure 3i. - Mass spectral scan of xenon isotopes 
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Figure 32. * Measured reaction product 
composition profiles for a methane -oxygen 
flame doped with SO 2 ♦ NaCI «H20. 
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Figure 33 • Comparison of measured and equilibrium calculated reaction 
products p'^oduced at 1413 K m a flowing gas mixture of \aCUa) $02igi 
♦ H^Otg- • O^ig' 
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